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EQUATION  OF  STATE  OF  STATIC  AND  DYNAMIC 
PROBLEMS  OF  THERMOPLASTICITY  UNDER 
COMPLEX  LOADING 

Yu.  G.  Korotkikh,  S.  M.  Belevich 

References  [l]-[8]  examine  the  coupling  equations  between 
stresses  and  strains  in  differential  form,  i.e.,  the  connection 
between  the  differentials  of  the  stress  and  strain  tensors.  In 
this  case  it  is  assumed  that  neutral  loading  (loading  tangentially 
to  the  surface  of  the  loading)  does  not  cause  plastic  deformation, 
but  the  hyperspherical  surface  of  loading  is  replaced  at  the 
point  of  loading  by  a hyperplane.  The  latter  is  exhibited  in 
the  fact  that  the  ccndition  that  the  new  surface  of  loading 
passes  through  the  terminus  of  vector  of  the  final  load  is  replaced 
by  the  condition  of  its  passage  through  the  terminus  of  the  vector 
which  is  the  projection  of  the  vector  of  final  load  on  the  direction 
of  the  gradient  to  the  surface  of  loading  at  the  point  of  loading. 

Experimental  determination  of  the  yield  surface  established 
the  specific  tolerance  to  plastic  strain  (•“0.002).  A decrease  in 
this  tolerance  means  that  the  fundamental  positions  of  the 
phenomenological  theory  of  plasticity  cease  to  be  valid.  This 
fact  leads  to  this  position  that,  on  the  one  hand  the  experimentally 
obtained  increases  in  the  stress  and  strain  tensors  cannot  bo 
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considered  as  differentials,  but  on  the  other  hand,  in  practical 
calculations  there  is  no  sense  to  determine  plastic  deformations 
to  within  less  than  tolerance,  i.e.,  also  it  is  not  possible  to 
consider  them  as  infinitesimal.  Furthermore,  in  practical  cal- 
culations, especially  during  the  solution  of  dynamic  problems, 
in  many  instances  it  is  difficult  to  regulate  the  size  of  the 
tensors  of  stress  and  strain  increments,  and  at  isolated  points 
they  can  reach  a significant  magnitude. 

This  work  attempts  to  deri/e  the  relationships  between 
increases  in  the  stress  and  strain  tensors  without  relying  on 
the  assumption  of  their  smallness,  and  gives  the  "physical"  | 
algorithm  of  practical  realization  of  these  relationships  when 
solving  problems  of  thermoplasticity. 


DERIVATION  OF  FUNDAMENTAL  PRINCIPLES 


During  derivation  of  fundamental  principles,  as  usual,  the 
stress-strain  tensors  and  their  increase  are  broken  down  into 
spherical  and  deviator  parts,  and  the  connection  is  established 
separately  for  each  part. 

The  complete  stress  tensor  e^j  is  broken  down  Into  elastic  e^ 
and  plastic  e^j  parts: 

£tj  ~ Sij  Cij.  | (l) 


At  a certain  point  in  time  let  the  stress  and  strain  state  of 
the  material  be  characterized  by  the  fields  of  the  stress-stain 


tensors  a 

_HOB  _ - 

cTjLj  ana 


£7  and  e™,  and  at  the  following  poi  . t in  time  by  tensors 

1h  OB  U 


'id 


Then  it  is  possible  to  write  the  relationships: 

.....  *• 

#-20<£v: 


(2) 


FTD-MT-24-2514-74 


2 


<■;?»,..  =»  2(0+  A 6) el/"",  | (J) 

where  Q*Q  (T,  o)  is  the  shear  modulus; 

AG  - change  in  the  modulus  due  to  change  in  the  temperature 

ai1 

T of  the  medium  and  the  hydrostatic  pressure  a=-^“-. 

By  deducting  from  equality  (2)  equality  (3)  and  carrying  out 
some  transformations,  we  will  obtain 


’ “ V •'**•••  ' *v  "A-'ij*  -u  ; ••*-».  *•  •»,>*  -•;*»»»?••  _ _ 


mob  y’ 


where  (see  Pig.  1) 


Ao'.  =*  2G  A <?'•?  + 

il  i ! 1 


AG‘‘  , 


(7  4-  AG 


- a. 

: •'/ 


HOB 


(4) 


r «°»y  — j' cr  4-  A a/. , ! 

//  <j  * ii  * ! 


:»<y 


Ae. 4-AeY ,+Ae?,  - strain  increment  tensor; 
ij  ij  ij 


>m,;j  || 


a*;;- 


= A-3,V  — Ac'.**  { for  case  a (); 

~ ( for  case  b j)  ; 
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Let  us  assume  that  the  tensor  of  increase  in  the  plastic  strain 
is  coaxial  to  a certain  tensor  Sjj , whose  sense  will  be  explained 
below? 


(5) 


By  utilizing  (*0  and  (5),  we  will  obtain? 


Contracting  in  (6): 


=»  2(G+A  G)  see  Pig.lt)  . ' 

a ....  w iu*  HxJi'  iii'iis . !u m&m? 


(6) 


we  will  obtain  'or  X? 


\ ; ■ 2(C+AG)S;,S*  ■ V 


(7) 


Let  us  explain  the  sense  of  tensor 


If  we  reject  the  previous  surface,  i.e.,  consider  that  the 
tensor  of  increase  in  plastic  strain  is  normal  to  the  previous 
surface  of  loading  at  the  point  of  its  intersection  with  vector 
Aa^,  then  by  S*^  it  is  necessary  to  understand  (see  Pig.  1). 

In  this  case  the  terminus  of  vector  Ac^  will  lie  at  the  point  of 
intersection  of  the  new  surface  with  the  line  passing  through 
point  c (terminus  of  vector  Ao^)  parallel  to  vector  S^j  (point  D). 
This  assumption  can  give  considerable  error  with  large  values  of 
increase  in  the  strain  tensor  Ae^CAo^j). 


In  this  case  it  is  expedient  to  accept  as  S 


ij 


the  vector: 


v.)  . . 

‘j 


m;  - v;,"'+  - o-®.)  %f\> 


where  6 and  6^  are  certain  constants: 

0<0<l,  OcSjd. 


when  6=61=0  the  vector  of  increase  in  the  plastic  strain  will  be 
normal  to  the  new  surface  of  loading  at  point  D.  With  0=0^=! 
we  will  obtain  the  first  case. 


The  connection  between  the  tensor  of  increase  in  permanent 
microstrain  and  the  tensor  .,f  increase  in  the  plastic  strain  is 
given  by  the  expression: 


where 


o-  A e?  4-  c o':?, 

c c j 1 1 1 j * 


(9) 


ft 


C = f ( | CU  I 3) 


is  determined  as  follows. 


K*-f  s,  T7re>iasp5tp,->  i ■ 


If  for  different  materials  we  use  experimental  data  to  graph 

IPijl 

as  a funotion  of  temperature  of  T,  where  p^j  are  coordinates 

c.f  the  center  of  the  yield  surfaoe  with  simple  elongation  at 
different  values  of  ep  and  temperature  T*0oC. 

Pi;j  are  coordinates  of  the  center  of  the  yield  surface  at 
certain  temperature  T (for  example,  see  Pig.  2),  then  accurate  to 
within  experimental  values  it  is  possible  to  consider  the  ratio 

iiHS 


M, as  independent  of  the  degree  of  plastic  strain.  This  graph 

j.fyl 

has  the  following  features: 


up  to  a certain  value 


.tot* 

T-T* 


*4 


const.*!,  then  it  sharply 


WM 

decreases  to  0 in  a narrow  temperature  range.  By  having  the 

M 


experimental  curve  JviU 

■i 


St; 


as  a function  of  T,  it  is  possible  to 


determine  c from  the  formula 


C = 


°r’  ^ when  A7’>0, 


( 9 *) 


where  csO  when  AT<0. 


The  first  and  second  terms  in  the  right  side  of  formula  (8) 
reflect  the  two  really  existing  coexisting  coexisting  processes: 
the  process  of  formation  of  microheterogeneities  with  plastic 
strain  and  the  process  of  recrystallization.  Temperature  T*  is 
the  temperature  of  recrystallization.  Here  we  do  not  consider 
the  fact  that  the  process  of  recrystallization  occurs  with  holding 
at  constant  temperature  (annealing). 


.•I 


!P?9S3gHS  W-W 
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" PHYSICAL"  ALGORITHM 


Cases  a and  b Pig.  1 are  not  distinguished  below.  Case  b 
differs  from  case  a in  that  In  case  a: 


'AC-iHi  = jA5;/!  - (A<n  ’ 

ana  in  case  b:  AC= 


The  iterative  process  of  determination  of  the  true  value  of 
increase  in  the  tensor  of  the  stress  deviator  {Aa^JsBD  can  be 
represented  by  the  following  algorithm: 


1)  In  the  beginning  we  take  0=8^1.  Through  the  terminus  of 
vector  {Ac.jJ}  (point  c)  we  draw  a straight  line  parallel  to 
(S*j)  up  to  the  intersection  with  the  previous  yield  surface, 

shifted  and  changed  because  of  AT,  a.,,,  e (point  D’)-  If  there  is| 
no  intersection,  in  any  manner  select  point  D1  (deflection  along 
normal  or  tangentially).  The  vector  of  stress  tensor  increment 
(AcjPj)  in  the  1st  approximation  is  defined  as  the  half-sum  of 
vectors  AD1  and  AC. 

Knowing  (Ac^)  in  the  1st  approximation,  we  compute  X 
according  to  formula  (7),  and  then  according  to  formulas  (5)  we 
compute  and  ^eij^l* 

We  compute  (S*^)  according  to  formula  (8)  at  some  values 
6 and  0^.  We  compute  the  new  value  of  the  radius  of  the  loading 
surface  (c^00)^^  as  a function  of  the  new  value  of  the  strengthening 
parameter,  calculated  in  terms  of  value  (AePj)^  according  to 
the  formula: 


(A’//>*  - 2G[Ae./  -,(A«f,),r+  — (’If),  f do 


8 


we  find  the  value  of  the  vector  of  the  stress  increment,  which 

We  test  to  see  whether  the 


corresponds  to  this  value  (Ae^,)^. 

» * 

terminus  of  vector  (a 


surface  with  center 


HOB 

id 


)2B<T^j  + (Ao^j ) g lies  on  the  new  loading 


=J?i  +(%)l  + cpff  I (11) 


radius  equal  to 


: V f i4r + 


This  condition  must  be 


fulfilled  accurate  to  a certain  e>0.  If  this  condition  is  not 
satisfied,  we  find  the  intersection  of  the  loading  surface  obtained 
in  this  approximation  with  the  straight  line  passing  through 
point  c parallel  to  (S*j)  (point  dhob1).  For  vector  C Aa ^ j ) 2 in 
the  second  approximation  we  select  the  vector  equal  to  the  half-sum 
of  vectors  (Aa^j ) and  vector  BDhob1  etc. 

The  procedure  for  the  correction  of  the  stress  tensor  proposed 
in  work  [9]  is  a particular  case  of  the  given  method  (p=0,  aT=const., 
e=61=0)  and  does  not  require  the  iterative  process.  If  there  is 
a yield  delay  effect,  the  given  algorithm  is  put  into  use  only 
after  fulfillment  of  the  conditions: 


V 


A ik 


t. 


> 1 , 


(11*) 


tr 

where  At  is  the  k-th  time  step,  and  t is  the  yield  delay  time 

°u-°T 

[ID],  which  corresponds  to  an  excess  of  i— at  this  step,  and 

dip  i 

K 

oy  is  the  value  of  stress  intensity  on  this  step.  In  this  case, 
after  fulfillment  of  the  conditions  (11),  it  is  possible  to  use 
different  laws  governing  change  in  the  dynamic  yield  point. 


The  above  relationships  between  the  strain  increment  tensor' 
and  the  stress  increment  tensor  can  be  presented  in  anisotropic 


9 


itzrm 


ggBggaaBMra«BH«B^^  •'  T 


form.  By  transforming  (5),  (7)  and  (4)*,  it  is  possible  to 
obtain: 


AeP 


'}£*•’■  rv  /-  *t  -r '■"■.■  ■ 


Aa^  . .*  V.  •• 

: I2CTaG) S-S-  7“  ” 2{G+4G) S'S*  “ ' 


• C* 


(12) 


(13) 


By  utilizing  then  the  equality  Ae^j=Ae^j+Ae^j , it  is  possible  to 
obtain  the  unknown  connection  between  the  increase  in  the  stress 
tensor  and  strain  tensor  in  anisotropic  form.  By  converting 
dependences  (12)  and  (13),  it  is  possible  to  obtain  the  expressions 
of  stress  tensor  increment  through  the  strain  tensor  increment. 


The  connection  between  the  spherical  stress-strain  tensor 
components  in  static  problems  can  be  presented  in  the  form: 


+ A(*n, 

3k  3k* 


(14) 


where  k=k(T)  -the  bulk  modulus,  a - the  coefficient  of  linear 
expansion,  Ae 


_Aei.1 
*“  • 


In  the  case  of  dynamic  problems  with  the  large  hydrodynamic 
pressures  which  occur  during  shock  wave  propagation,  it  is 
necessary  to  utilize  the  equations  of  state  of  solid  bodies  at 
high  pressures  and  temperatures  in  order  to  know  the  connection 
between  spherical  components. 

As  shown  in  work  [11],  a large  role  in  these  equations 
belongs  to  the  dependence  of  pressure  on  internal  energy. 
Subsequently,  the  connection  between  the  spherical  stress-strain 
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tensor  components  is  established  on  the  basis  of  equations  of 
state  of  Mie-Grueneisen  [12-133 


P 


(15) 


where  Pn(v)  - the  potential  part  of  the  pressure,  caused  by 
compression  of  crystal  lattice  at  0°  Kelvin.  This  dependence 
can  be  obtained  from  the  experimental  data  [lM]  in  the  form: 


p--S* 

i 


(16) 


where  ai  - certain  experimentally  determined  coefficients,  and 
p - density  of  element  (p/Pq)=( Vqv) , ^(E-Ev)=PT  - thermal  part 
of  pressure. 


Y=y(v)  - Grueneisen's  coefficient.  It  is  given  in  the  form: 


where  A,  B,  C are  experj mentally  determined  coefficients. 


If  the  shock  adiabat  is  knov/n  for  just  one  material  with 
normal  and  low  density  (porous),  then  y(v)  can  be  determined  by 
the  formula: 


2v(Pq  - j;  : 

ll  '■  - — 

PqVqq  — Pvu  ’ 


(18) 


or  calculated  from  certain  approximation  formulas  given  in  [16], 


E - specific  internal  energy. 


! 

11  : 


•*•,  -.*.-*  . — . , » . •■;  ■*;  t~:^  ,+  t.,  •^vy*-  • •*  * \‘;*  ?**;^V,V/:'\  ***'■:  !?%v  ''?*!"  •,*'‘!'  v " ■•■^‘x  -. 

\e  = EJiV)  + Et  = — j p(v)  dv  + Aaf“  + j cvdT,+En , :i9) 


where  Azc,ab  - energy  of  deviator  strain,  and 


'v  - specific  heat  of  lattice  at  constant  pressure. 


To  equation  (15)  it  is  necessary  to  add  a certain  thermodynamic 
identity,  which  characterizes  the  strain  condition.  This  can  be 
either  condition  on  the  discontinuity  - the  Rankine-Hugoniot 
equation 


v ';\|  * “■*  v-v'?^  y;  vtr  • r NJ  * 

E-E»--~  ~ ( P + p0)  (®c-»)  + A 

2*  . -.ACVjCi.?'. 


(20) 


or  the  equation  of  the  adiabatic  curve 


• *?s  rs 

■ dS  = o;; 


(21) 


where  dS  is  the  increase  in  the  entropy,  or  any  other  thermodynamic 
relationship,  which  expresses  the  condition  of  conservation  of 
energy.  The  use  of  the  indicated  relationships  makes  it  possible 
to  consider  the  fact  that  the  shock  loading  of  a material  passes 
along  the  shock  adiabat,  and  unloading  passes  along  the  isentrope. 

Utilizing  equations  (15),  (20)  or  (21),  it  is  possible  to 
calculate  the  change  in  temperature  at  a particular  point  of  the 
body  [12];  temperature  is  determined  from  the  integral  relationship. 


(22) 
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where  FQ  is  initial  thermal  energy:  £*)  — j C p l % Jand  cy  is 
determine^  from  formula  [12]:  .0..  - “ ' 


Cv 


35/  T | 

e*‘T  - 1 1 

...* ** 


Is 


(23) 


R - the  universal  gas  constant, 
y - molecular  weight, 

D - Debye  function 
6 - Debye  temperature. 


Since  at  T=0,  cy  reaches  961  of  its  maximum  value,  when  T|.e  cy 
is  considered  constant,  and  Debye  function  D will  take  the  form 


Energy  sources  of  neutron  irradiation  type  are  considered  through 
the  increase  in  internal  energy  of  the  element.  After  determining 
the  connection  between  the  spherical  and  deviator  stress-strain 
tensor  components,  we  obtain  the  complete  relationships  between 
stresses  and  strains,  which  determine  the  behavior  of  the  medium 
with  nonisothermal  elasto-plastic  deformations. 
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